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We investigate the distribution of low-lying eigenmodes of hermitian Wilson-Dirac operator, Hw, 
with the gauge action whose form is designed to avoid topology change, as well as with the 
standard plaquette action. On the quenched lattices, the former gauge action exhibits less density 
of low-lying eigenmodes of compared to the latter at the same lattice spacing. We also show 
preliminary results for the dynamical simulation with two flavors of overlap fermions. 
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1. Introduction 

Recent progress of the chiral fermions, both in theoretical understanding and in technical de- 
velopment of numerical algorithms, is making its dynamical fermion simulation more feasible with 
the next generation machines. Wide range of applications from flavor physics to hadron physics is 
expected with the dynamical fermions having exact chiral symmetry. 

The Neuberger's overlap-Dirac operator with mass m is written as [|l]] 

An-H = (l - ^) ^ov + m, Dov = ^[1 + 75sign(//iv)], (1.1) 
where a = a/{\+s). For the kernel Hy^ we use the Wilson-Dirac operator Dyi/ as 

Hw = y5[Dw-{l+s)]. (1.2) 

The parameter s is chosen in the range 1^1 < 1. The locality of the overlap operator is proven when 
the bulk of eigenvalues of Hy/ is bounded in a non-zero region [Q]. This condition is satisfied 
when the gauge configuration is sufficiently smooth. In practice, however, it is known that there 
is non-zero density of near-zero eigenvalues when the gauge configurations are generated with the 
conventional Wilson gauge action |^. Although the problem of locality is avoided if the near-zero 
modes are localized [Q, ||, ^, it is more desirable if the gauge action does not generate such near- 
zero modes from the beginning. The use of improved actions (LUscher-Weisz, Iwasaki, DBW2) is 
a popular choice, with which the density of the near-zero modes is reduced, but in this work we 
consider an alternative solution. 

Since the zero mode of Hw appears when the topological charge of the gauge field changes 
its value, we expect that the appearance of the near-zero mode is suppressed if we use a gauge 
action for which the topology change is suppressed or even forbidden. In this work, we employ 
the Liischer's topology conserving gauge action [Q] and investigate how much it suppresses the 
near-zero modes. 

The suppression of topology change and of the appearance of the near-zero modes is desir- 
able also in a practical point of view. (1) The cost of approximating the sign function in ( |1.1[ ) is 
reduced; (2) In the dynamical fermion simulation, one can avoid the "reflection-refraction" process 
[^], which is very costly; (3) Gauge configurations with a fixed topological charge is efficiently 
accumulated, which is useful for the study of the £-regime of chiral perturbation theory. 

2. Topology conserving gauge action 

The topology conserving gauge action can be constructed by enforcing the admissibility con- 
dition 

\\-P^v{x)\<e, yx,n,v, (2.1) 

where P^vix) is the plaquette. £ is a parameter to control the smoothness of the gauge field. For 
the SU(3) gauge theory, the locality of the overlap-Dirac operator is guaranteed when e < 1/20.49 
[^]. In practice this condition is too strong to realize the lattice spacing ~ 0.1 fm and we have to 
relax it to £ ~ 0(1). Therefore, the question is whether the good property is kept enough for these 
large value of £. 
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Figure 1: Low-lying eigenvalue distribution on quenched 16'* lattices for 1/e = 0, 2/3, and 1 (a ' ~ 2.5 
GeV). 



An example of the gauge action to enforce (|2. 1[) is 



Sg 



1 



Sp{x) 



^^,{l-Sp{x)/eY 



(2.2) 



where Sp{x) = 1 — ReTr/'^v(-^)/A'f ■ In this work, we set a 
used by two of us for the study of the massive Schwinger model [10]. In the four-dimensional 



if admissible 
otherwise 

: 1. Previously, this gauge action was 



QCD, Bietenholz et al. are also studying its property [11, 12]. 

In the following, we investigate the eigenvalue distribution of H-w on the gauge configurations 
generated with this action. The topology change and scaling property of this action is presented in 
these proceedings [p^. 

3. Quenched results 

First we examine the low-lying eigenvalue distribution of H-\/v for the topology conserving 
gauge action in the quenched approximation. The configuration is generated with the Hybrid Monte 
Carlo algorithm with 8t = 0.01 and a length of trajectory 0.2. The parameter £ is set to l/£ = 0, 
2/3, and 1. (l/£ =0 corresponds to the standard plaquette action.) The admissibility condition, 
1 — 5'p(.x)/£ > 0, is monitored during the molecular dynamics evolutions, and no violation has 
been found for 1 /e =2/3 and 1 . For finite l/£, while the topology change is not exactly prohibited, 
it is indeed much suppressed compared to the plaquette action 1 /£ 



0[13] 



Figure H shows ten lowest-lying eigenvalues of \Hy/ \ on a 16^ lattice for the three values of £. 
The inverse lattice spacing determined with the Sommer scale is oT^ ~ 2.5 GeV for each lattice. 
In the plot, A, B, and C correspond to the combinations (j3,l/£) = (6.13,0), (2.70,2/3), and 
(1.42,1), respectively. As expected, we find that the eigenvalue distribution shifts upward as l/£ 
increases. 

Figure ^ shows the results on 20^* lattices. In this figure, the low-lying eigenvalue distribution 
is compared at five sets of parameters: (j8, 1 /e) = (6.0,0), (2.55,2/3), (1.3,1), (2.7,2/3), and (1.42,1). 
The first three correspond to a^' ~ 2. 1 GeV, while the last two are 2.6 GeV and 2.5 GeV. The same 
tendency as in Figure |l] is found. 
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Figure 2: Low-lying eigenvalue distribution of IH-w on quenched 20^ lattices. 
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Table 1: Parameters of our dynamical simulation runs. 



We conclude that in the quenched calculation the density of the near-zero eigenvalues is 
smaller for the topology conserving gauge action than that of the plaquette action. This trend 
is enhanced as 1 /£ increases. 

4. Dynamical results 

In this section, we show the results of an exploratory study of the low-lying eigenvalue dis- 
tribution of Hy^ in the dynamical fermion simulation with two flavors of the overlap fermion. Our 
algorithm of the Hybrid Monte Carl update follows that of ||8|]. Namely, we adopt the Zolotarev 
rational approximation [ p4| | with degree Np = 20 to approximate the sign function. The low-lying 
eigenmodes of Hw is calculated using the implicitly restarted Lanczos method. When the low- 
est eigenmode changes its sign, the value of the pseudofermion action changes discontinuously. 
Then the reflection or refraction occurs depending on the energy of the momentum of the mode 
parpendicular to the surface corresponding to that eigenmode. The occurrence of the reflection and 
refraction is monitored during the HMC update. 

Numerical simulation is performed on 4^ lattices, with overlap fermions at m = 0.2 and s = 0.6. 
Configurations are generated with 5t = 0.02 and an unit length of trajectory. The values of jS and 
l/£ are listed in Table [I| Since the lattice spacing is unknown at this stage and the lattice volume 
is so small, we cannot draw any definite conclusion, but we would report the status of our study. 

The average number of occurrence of reflection is listed in Table and the refraction is used to 
trace the change of topological charge, as shown below. If the topology conserving action properly 
works, appearance of the near-zero mode should be suppressed, and hence the reflection/refraction 
should also be suppressed. Our result does not show clear tendency supporting this expectation. 
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Figure 3: History of plaquette (top), topological charge (middle), and the low-lying eigenvalues (bottom). 
Results for the plaquette gauge action (1/e = 0) at j3 =5.40 (left panel) and j3 = 5.50(right). 



The simulation time history of the plaquette, the topological charge, and the low-lying eigen- 
values of //ly is displayed in Figures ^ and ^ for 1/e = and 1, respectively. The topological 
charge is observed by monitoring zero-crossing of the lowest eigenvalue of Hw. No clear differ- 
ence is found for the topology change between the plaquette and the topology conserving gauge 
actions. For the distribution of the low-lying eigenmodes, the topology conserving gauge action 
exhibit smaller density than the plaquette action. 

To draw definite conclusions, we of course need to compare these actions on fixed lattice 
spacings on larger lattices. The present runs are probably done on rather coarse lattices. Whether 
the topology conserving action is practically feasible or not will be examined in more detail. 

This work is partly supported by the Large Scale Simulation Program No. 136 (FY2005) of 
High Energy Accelerator Research Organization (KEK). H.M. is supported by Grant-in-Aid of the 
Ministry of Education (No. 16740156). 
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Figure 4: Same as Figure 1 but for 1/e = 1: j3 = 0.70 (left panel) and /3 = 0.80 (right). 
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